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An easily accessible method is presented that permits to calculate spectra involving atomic mul- 
tiplets relevant to X-ray Absorption Spectroscopy (XAS) and Resonant Inelastic X-ray Scattering 
(RIXS) experiments. We present specific examples and compare the calculated spectra with avail- 
able experimental data. 
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I. INTRODUCTION 

Since the early days of quantum mechanics and atomic 
spectroscopy, it has been realized that electron-electron 
interactions and spin-orbit coupling were at the origin of 
the splitting of the electronic shells into multiplet levelsX 
As they originate from a many-body problem, these levels 
are relatively straightforward to compute and understand 
analytically as long as the number of involved particles 
remains small. Relatively rapidly one is forced to require 
computer assistance. To this effect, Cowan2, implemented 
a systematic computer code in order to solve multiplet 
structures arising in atomic spectra. 

Core spectroscopics have gained interest in solid state 
physics in the eighties as it has been shown that they can 
provide deep insights about the electronic structure of the 
materials. As the processes involve one or more open- 
shell they exhibit signatures of multiplet structures ,^^1^ 
these structures being sensitive to the presence of local- 
ized moments and hybridization with the ligands. Later 
on, using a multiplet-based code, a systematic study of 
XAS in 3d transition-metal compounds within a cubic 
crystal field environment as been performed by de Groot 
et al.^ showing good agreement with experiments. In- 
deed, soft x-ray XAS^ and also RIXS^ are linked to lo- 
cal physics: the absorption of a photon and creation of 
a localized core-hole opens up a shell and therefore the 
multiplet structure becomes apparent in the spectra. 

With the development of sources and progress in op- 
tics. X-ray spectroscopies such as XAS and RIXS, espe- 
cially in the soft X-ray regime;^ became tools of choice 
to investigate transition correlated materials such as, for 
instance, cupratesiSiii or vanadates ji^J^ii^ the current 
resolving power allows to investigate multiparticle exci- 
tations like it was predicted in different theoretical con- 
texts for inelastic X-ray scattering] ^^i^^i^^i^^ From here 
we see that it becomes crucial, while interpreting the ex- 
perimental data, to have a systematic, user-friendly and 
transparent way of computing the multiplet spectra in 
order to disentangle in the experiment the information 
arising from single-particle excitations from the informa- 
tion relevant to collective excitations. 

As noticed earlier, atomic multiplets code already ex- 
ist in the X-ray spectroscopy communityr^iSSii^ they are 
based on Cowan's pionneering work and make use of sym- 
metries in a very educated way. However, for a gen- 
eral audience and especially for newcomers to the field 



of X-ray spectroscopy, we felt that with nowadays com- 
puting power we had the opportunity of writing a less 
expert-oriented code, in that sense that the diagonal- 
ization of the Hamiltonian can be done brute force and 
the local crystal field can be implemented in a straight- 
forward manner. Atomic-multiplet calculations for ar- 
bitrary symmetries have already been reported^ how- 
ever the implimentation of the crystal-field in the latter 
approach still requires the introduction of free parame- 
ters. The aim of this paper is both to present our ap- 
proach which combines the Density Functional Theory 
(DFT) -based solution form the radial Dirac equation 
to a full diagonalization procedure in order to compute 
the multiplet structures as well as the associated XAS 
and RIXS spectra by treating all interactions (electron- 
electron, spin-orbit coupling and crystal-field) on the 
same footing. 

The paper is organized as follows: in a first part we 
remind the reader about atomic multiplets as described 
extensively in Refs.[lll^l7ll23j. we discuss how to imple- 
ment the crystal field in our calculation. From these ba- 
sic steps we are in a position to compute explicitely XAS 
and RIXS spectra ; two sections are devoted to these 
issues where we also discuss the implementation of the 
optical selection rules in the dipolar approximation. 



II. ATOMIC MULTIPLETS: OVERVIEW 
A. Electronic interactions and Hilbert space 

By solving the one-particle Dirac equation we know 
that the electronic states are determined by quantum 
numbers n and j. In a non-relativistic treatment we 
can define an orbital quantum number I and shells for 
the orbitals (like s, p, d, / for respectively I = 0, 1,2,3) 
such that these shells have a {21 + l)-fold orbital degen- 
eracy times the Kramers degeneracy for the spin degree 
of freedom. For instance, we see that (2x2 + l)x2 = 10 
electrons can fit in a c?-shell. 

If we consider the case of a non-hydrogenic atom, with 
more than one electron, the electronic eigenstates and 
eigenvalues are not simply labeled by the quantum num- 
bers n and I because the electrons interact electrostati- 
cally with each-other (and also via spin-orbit coupling). 
As long as a shell is empty or fully occupied it remains 
symmetric and the orbitals of a given shell remain de- 
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generate. However, for an open shell electrostatics and 
spin-orbit interactions will split the shell and form a mul- 
tiplet structure. 

Hence, the Hamiltonian describing the system of elec- 
trons is given by : 



(1) 



where the first term is a single-particle diagonal operator 
containing the kinetic and automatically the spin-orbit 
energy since we use solutions from the Dirac equation. 
The second term is simply the sum over all pairs of elec- 
trons for electrostatics interactions. This last term being 
off-diagonal we see that we end up with a diagonalization 
problem in order to deduce the niultiplet structure. 

In the present paper we are aiming at describing soft 
X-ray absorption and emission processes involving two 
shells. Hence, we can define the Hilbert space associ- 
ated to this problem by only considering the electrons 
involved in these two shells. For a shell determined by 
the quantum numbers (n, /) and containing k electrons, 
the dimension of the Hilbert space is : 



M{l,k) 



[2(2/ + 1)]! 



[2(2/ + 1) - k]\ k\ 



(2) 



For instance 2 electrons in a p-sliell gives A/'(2, 2) = 15 
states, 2 electrons in a d-shell gives A/'(3, 2) = 45 states. 
If more than one shell is opened, each shell being indepen- 
dent, the overall size is given by the product A/" = Hi 

While the case of transition metals can be handled 
on standard desktop computers, it is worth mentioning 
the case of lanthanides. For instance, the largest Hilbert 
space is given for the case of Gd Af^bd^ — *■ AphcP 
and has a dimension 135135 for its final-space, it ne- 
cessitates a too large amount of memory. Fortunately, 
most of the Hilbert spaces involved for soft X-ray scat- 
tering experiments are of the order of Sm at M4^5-edg6^ 
Ad}^Af^ 4(i^4/^ which are of respective dimensions 
3003 and 34320 (^ 9GB), such sizes being available on 
modern computers. 

Care has to be taken while generating the Hilbert 
space ; the space is obtained by taking the direct product 
of the one-particle states in the Fock space, yet, electrons 
are fermions, therefore they obey Pauli exclusion princi- 
ple and they anticommute. 



B. Lattice and Crystal Field 

So far, Eq.([T]) did involve the matrix elements for a 
multiplet formed by the electron-electron interaction and 
the spin-orbit coupling of a single isolated atom, but we 
have to keep in mind that the aim of this paper is to 
simulate multiplet structures of ions in solids. Therefore, 
we consider in this section an ion embedded in a clus- 
ter; the considered ion feels the electrostatic crystal field 
potential created by the other ions forming the cluster. 



TABLE I: Positions of the ligands for a 4% uniaxially dis- 
torted perovskite. The position are in units of the typical 
distance d =1.96A, the formal charges are in units of e. 
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charge 


-1.00 


0.00 


0.00 


-2.00 


1.00 


0.00 


0.00 


-2.00 


0.00 


-1.00 


0.00 


-2.00 


0.00 


1.00 


0.00 


-2.00 


0.00 


0.00 


-0.96 


-2.00 


0.00 


0.00 


0.96 


-2.00 



The crystal field splitting in simple geometries (cubic 
or its subgroups) can be deduced from symmetry argu- 
ments. In our approach we choose to implement the 
crystal field in the fuU-diagonalization procedure. For 
the implementation of the crystal field we assume a point 
charge model : the ion with core excitation is surrounded 
by point charge ligands. 

The advantage of solving the problem of a cluster con- 
taining the ion and the neighboring ligands lies in two 
facts: i) there is no need to start from a cubic symmetry 
but the position of each anion can directly be entered 
in the program ; ii) the crystal-field strength is not a 
free-parameter but the code allows to be predictive. 

The potential created by the ligands of charge qi at a 
distance Ri on the orbital r is given by the electrostatics : 



N 



(3) 



assuming that Ri > r the previous expression can be 
expanded on the spherical harmonics basis, as we will 
see in the appendix (Eq lAip . 

The introduction of a crystal field in the problem 
explicitely breaks the spherical symmetry of the prob- 
lem, therefore we have to choose a local cartesian basis 
(ex,ey,ez;). A typical distance d (in A) is introduced 
defining the unity of the coordinates system (x,y, z), 
such that the surrounding ions are represented by their 
positions in this basis and the formal charges that they 
carry in units of e. 

A straightforward example is a d^ ion (without spin- 
orbit coupling) in the perovskite structure. Consider- 
ing the octahedral environment only, for the undistorted 
case, two levels should be observed while, with a slight 
distortion of the octahedron, the degeneracy of the 
level is partially lifted and the two Cg levels split-up as 
well. The schematic shown in Fig[T] defines the (x, y, z) 
basis chosen to define the local crystal field such that in 
units of d, the 6 ligands are located at (±1, 0, 0), (0, ±1, 0) 
and (0, 0, ±dz/d). The positions and charges of the ions 
corresponding to this case is summarized in Table |T] and 
reflects the input parameters of the code. 

The case of a distortion in a perovskite structure of 
the type of LaTiOa is displayed in Fig. [21 It is relatively 
easy to understand qualitatively what is happening in 
Fig. [21 for the cubic symmetry, the splitting is observed 



FIG. 1: (Color online) Ion in a distorted octahedral environ- 
ment; the ion with core excitation is in grey while the red 
spheres represent the ligands, in green the incoming photon 
with its associated polarization vector e. 
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FIG. 3: (Color online) Plot of the radial part of the d-wave 
function for Ti as a function of r in atomic units : the tail of 
the wave function shrinks with the ionization. 



Compared to a summation over the whole lattice, the 
inclusion of the crystal field by only considering the sur- 
rounding environment of an ion remains relevant for the 
present investigation. Indeed, we can include many ions 
-and not only the closest ones- such that the symme- 
try of the crystal field is preserved. Once the symme- 
try independent parameters of the crystal field are ex- 
hausted, addition of further more neighbor shells has a 
minor quantitative effect. Screening of the bare crystal 
field may bee taken into account by a scaling factor. 



Effect of Ionization 



FIG. 2; (Color online) Partial lifting of the degeneracies as 
an effect of the distortion, the orbital degeneracy is indicated 
in parenthesis. Inset: Undistorted and distorted perovskite 
structure in blue Ti''"'', in red and in green La^^ ions. 

because the eg orbitals {dx2-y2, d^^'^^r^) have their lobes 
pointing towards the 0^~, therefore they are energeti- 
cally disfavored compared to t2g orbitals. When the dis- 
tortion takes place, the two apical oxygens are moved 
closer to d^z'^_j.2 orbital on one hand, and to dxz, dyz or- 
bitals on the other hand, leading to the splitting of both 
sub-shells. 

It may be of interest to notice that traditionally, the 
crystal field splitting between t2g and eg levels is referred 
to as lODq, where the parameter Dg is defined by con- 
vention : 

6Dg = JriegWi^ieg) dr 
-4Dg = Jr{t2g)Vyj{t2g) dr 

such that the baryccntcr of the splitting lies at 0. Dq is 
thus taken as a measure of the crystal field strength and 
often adapted to fit experiments. 



The ionization has a major impact on the radial part 
of the wave-functions. As a consequence of taking one 
or many electrons out of an atom, the remaining elec- 
trons feel a less screened nuclear potential and the wave- 
functions contract around the nucleus, this effect can eas- 
ily been seen in Fig. [31 

To be quantitative, for this specific case, the expec- 
tation value (?'), has a variation of about 30% between 
the neutral titanium ion and Ti^+ ion. At first glance, 
this does not seem to be dramatic but this quantity 
enters in powers of / in the evaluation of the crystal 
field, and we see that taking formal charges too seri- 
ously can have a strong impact leading to underestimat- 
ing the crystal field strength. To avoid this situation 
we consider systematically neutral ions while computing 
the radial wave-functions and allow for a scaling of the 
crystal field strength if necessary. Indeed, self-consistent 
band-structure calculations suggest that neutral atoms 
represent the charge density around an atom more re- 
alistically than ionic. We checked on different transition 
metal oxides (at the L-cdge) that this approximation pre- 
dicts reasonably well crystal-field strengths and docs not 
necessitate the introduction of other fitting parameters. 
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D. Spin-orbit coupling 

Previous atomic multiplet calculations, even recent 
onesr^ often considered solutions of the nonrelativistic 
Schrodinger equation, meaning that the spin-orbit cou- 
pling had to be introduced by hand in the Hamiltonian 
matrix. Since we are willing to treat all interactions on 
the same footing, the most straightforward approach is 
to consider the solutions of the Dirac equation. Hence, 
the spin-orbit coupling manifests itself in two ways: in 
Eq.([T]) the values depend on the total angular mo- 
mentum J as well as the radial wave-functions do. From 
our experience and as we will show in the next section, 
for X-ray absorption spectroscopies, this treatment of the 
spin-orbit coupling allows a relatively good first-principle 
prediction of the L2-L3 edge splitting. Yet, we allow for 
an arbitrary scaling of this coupling which enables the 
user to have a better match with experimental data if 
needed. 



E. Polarization 

As noticed earlier, the introduction of the crystal field 
breaks the spherical symmetry of the pure atomic case. 
The description of X-ray scattering experiments necessi- 
tate to discuss in this subsection the implications of the 
crystal field on polarized photon experiments. The in- 
coming photons are defined by both their energy (huj) 
and polarization e. the polarization property plays a 
central role in the scattering process: depending on 
the polarization a transition is dipole-allowed or dipole- 
forbidden as we will show in the next section. Thus, it 
is crucial to implement the polarization in the problem. 
The most natural way of introducing the polarization is 
by referring to its projection [cx, Cy, Cz) on the same lo- 
cal basis (x,y, z) which has been defined for the posi- 
tion of the surrounding ions. For the specific example of 
FiglU the polarization associated to the incoming photon 
is e = (0,0,1). 

The particular case of unpolarized photons, or exper- 
iments done on powder samples, we can always recover 
this limit by doing an incoherent superposition of the 
different polarizations. 



III. X-RAY ABSORPTION SPECTROSCOPY 
(XAS): 

A. Description 

X-ray absorption spectra are obtained by shining an 
X-ray photon on a material, the photon is absorbed and 
gives rise to a transition of a core-electron into an excited 
level. Thus a core-hole is formed and therefore at least 
one shell is opened resulting the formation and investi- 
gation of the multiplet structure. 



The absorption spectra can be evaluated with the 
Fermi Golden rule : 

I{u;) = mOHo)\'S{c^ + Eo- E,) ; (5) 

i 

where (|^/'o)i ^-o) refers to the ground-state and its eigen- 
value, {\tpi), Ei) to the final state and final energy, 
whereas uj is the energy of the absorbed photon. The 
operator O represents the transition operator which is 
treated in the dipolar approximation. 

Equation ([5]) has its limitations both from an intrinsic 
and experimental point of view : indeed the core-hole 
has a finite intrinsic lifetime and the experiment has of 
course a finite resolution. These two facts lead to spec- 
tral broadening which we can mimic by a Lorcntzian and 
Eq.® becomes : 

'v^^ (6) 

i ^ ^ 

In practice, the broadening is partly due to the exper- 
imental resolution but it should be mentioned that the 
broadening term F may not remain a constant but can be 
a function of uj and ■(/'i '■ this reflects hybridizations and 
vibrational effects which depend on the orbitals involved 
in the states as well as electron-electron scattering. 

To keep a straightforward approach to our calculation 
we will assume a finite but constant broadening which 
can be tuned to match the experimental limitations. The 
most striking consequence of this assumption lies in the 
fact that the relative intensities between the peaks in our 
simulations may be slightly different from what would 
experimentally be observed. 

B. Dipolar Approximation 

The incoming photon is defined by its polarization e 
and its wavevector k. If we assume that the vector po- 
tential can be expanded in plane- waves, keeping in mind 
that k • r <C 1 we can write the matrix elements we need 
to compute for the transition like : 

(V'.le • pe^'^'lVo) - me • p (1 + ik • r - ...) |^o) (7) 

In the electric dipole approximation we replace the ex- 
pansion by 1, such that the operator is given by : 

6 = e • p = e • [r, cx e • r (8) 

Using an expansion on the y/", we finally obtain : 

cx r (eiYii + £0^1° + ^-lY^^) (9) 

where the coefficients Ci represent the projection of the 
polarization vectors on the Y™ basis. 

At the dipolar approximation level, there are some 
selection rules for the transition which involve different 
quantum numbers : 

a; = ±1 ; As = ;AJ = 0,±1; ... (10) 
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FIG. 4: (Color online) In red and blue 2p'^3d^ multiplet en- 
ergy levels. The length of the bars is proportional to the 
degeneracy (in red J — 1). In black XAS L-edge spectrum 
for a Ti*""" free-ion : we see that only the final states with a 
AJ = 0, ±1 with respect to the ground state contribute to 
the spectrum. 

It should be mentioned here that these rules are auto- 
matically included in our code through Clebsch-Gordan 
and Gaunt coefficients when we compute the overlap. 

As an example, in Figd] we take the standard case of 
an L-edge absorption for a Ti^+ ion : 2p^3dP 2p^3d^. 
The effect of the dipolar selection rules is clearly seen : 
since the ground-state \tPq) has a total momentum J = 0, 
only the three states with J = 1 in the multiplet of 
contribute to the formation of peaks in the XAS spec- 
trum. The inclusion of a crystal field would change the 
character of the states in the multiplet, allowing there- 
fore for new dipole transitions and the appearance of new 
peaks. 

C. Comparison to experiments : ATiOa 

ATiOs (A=Sr,La) are perovskite structure materials: 
the Ti ions are in a cubic environment. Since Sr and 
La carry a formal charge 2+ and 3-1- respectively, at 
the titanium L-edge the optical transitions are given by 
2/3d° 2p^3d^, and 2^63^1 ^ 2p^3(f respectively. 
The corresponding absorption spectra are displayed in 
Figs El and [H] 

It is relatively straightforward to understand the 
2/3d° 2p^i(f XAS spectrum of FigE) the two 
groups consisting of two large peaks correspond to the 
L3-L2 splitting while the peaks themselves represent t2g 
(~ 458eV" ; ~ 463. SeF) and Cg (~ 460ey ; ~ 465. SeV) 
levels. The two small satellites close to ~ 457eV are due 
to the fact that the final states are in fact a bit more com- 
plicated than this simplistic view : to the t2g-eg splitting 
occurring for the 3d} configuration we should also con- 
sider the effect of the 2p^ shell, which will give rise to 



FIG. 5: XAS spectrum for a purely 2/3d"-case: SrTiOa at 
the titanium L-edge. Upper panel: experimental result ob- 
tained by J. Schlappa et al. Ref. [S^]. Lower panel: multiplet 
calculation 
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FIG. 6: XAS spectrum for a 2p''3d^-case: LaTiOa at the tita- 
nium L-edge. Upper panel: experimental result obtained by 
T. Higuchi et al. Ref. [231. Lower panel: multiplet calculation 



these features. 

For the 2p^3d^ 2p^3d^ transition we still can observe 
a L3-L2 splitting on Fig[6l as well as a t2g-eg splitting al- 
though it is less clear than in the previous case since here, 
the effect of electron-electron interaction tends to give 
rise to a more complex multiplet structure. Furthermore, 
LaTiOs is known to be a Mott-Hubbard insulator^ for 
which the electronic correlations play a more important 
role than in SrTiOs, therefore the present fully ionic ap- 
proach partly fails to reproduce the XAS spectrum in a 
truly quantitative way. 

Although our present approach does not take charge 
fluctuations explicitly into account we still can mimic 
their impact on an absorption spectrum whenever the 
charge fluctuation is relatively weak. To this effect and as 




FIG. 7: XAS spectrum for a doped-case. Upper panel: exper- 
imental result obtained by T. Higuchi et al. Ref . [27| . Lower 
panel: multiplet calculation arising from the superposition 
+0.1/^1. 



example, we consider in this paragraph a shghtly doped 
case for these titanates compounds: Sro.gLao.iTiOa. The 
wave-function for the ground-state will then consist in a 
linear superposition of the two previous one : |^o} = 
0.913^") -I- O.ljSd^). The result is compared to the ex- 
perimental data taken by T. Higuchi et al. in Ref. [27| . 
Although the writing of the wave-function as a super- 
position is formally correct, the evaluation of the XAS 
spectrum from a superposition remains an approxima- 
tion, neglecting interference terms. The expression of 
Eq.(l5]) would contain interference terms which vanish in 
the present approach. The approximation gives a base 
to describe the experiment, by comparing to the un- 
doped case of FigO we can see that the effect of the 
doping largely contributes to the formation of satellites 
peaks around the main features which were previously 
addressed to Cg and t2g levels, this effect leads to a change 
in the relative weight and intensity of each peak. Thus a 
direct and systematic comparison of the computed XAS 
spectra for different doping levels can give insights on the 
respective role of the Cg and t2g bands that are formed 
and contribute to experimental XAS data. 



IV. RESONANT INELASTIC X-RAY 
SCATTERING (RIXS) 

Resonant Inelastic X-ray Scattering (RIXS) is a 
photon-in photon-out spectroscopy with an incoming 
photon energy hiOin close to an absorption edge, such 
that the denominator of the second order-term in a Fermi 
golden-rule is small and therefore completely dominant. 
Being of second order nature, we see immediately that 
the RIXS process involves not only an initial j^'o) and a 
final \Tpf) state but also an intermediate state \ipi) with a 
finite core-hole lifetime F,- which can be taken as a con- 



FIG. 8: (Color online) Upper panel: RIXS spectrum on 
MnO taken at the Mn L-edge, data by G. Ghiringhelli et al. 
Ref. [2^. Middle panel: spectrum obtained with the multi- 
plet approach for unpolarized outgoing photons. Lower panel: 
spectra obtained with the multiplet code. The different colors 
are assigned to different polarizations of the outgoing photon: 
in-plane and in vertical polarizations. 



stant. 

In a first time, starting from the ground-state of energy 
Ei, the intermediate state, at energy is accessed via 
an optical excitation and hence depend on the polariza- 
tion the incoming light ej„. The second part of the pro- 
cess consists in the de-excitation from the intermediate 
state to a final state of energy Ef, in emitting a photon 
of energy hujf. This process results in the writting of the 
Kramers-Heisenberg formula : 



/ 



^ Ei ~ Eo ~ hcut 



6 {Ef - Eo- nn) 

(11) 

where the same optical selection rules as in the previ- 
ous paragraph will apply. In the later expression we 
see that the intensity also depends on the energy-loss 
hfl = h{ujQ — LOf). The experimental RIXS spectra also 
contain broad fiuorescence lines which mainly come from 
the emission from intermediate states with finite lifetimes 
which are not taken into account in the multiplet cal- 
culation. It should be mentioned that, apart from the 
fluorescence, the broadening of the ^-function, in the en- 
ergy loss direction, depends on the lifetime of the final 
state IV'/) and on the experimental resolution, whereas 
the broadening in the incoming photon energy direction 
is mostly governed by the lifetime of the core-hole inter- 
mediate state. 

In order to have a specific example we compare with 
the Mn L-edge data by Ghiringhelh et al. for MnO.^^ 
From a structural point of view, MnO has the same crys- 
tal structure as NaCl, on the charge distribution side, 
although MnO is not purely ionic, our code describes a 
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static Mn^+0^~ picture, however we can adapt and scale 
the crystal field by changing the formal charge of the lig- 
ands. Hence, at the L-edge, we have transitions of the 
type: 2^63^5 ^ 2^53^6 ^ (2/3^^)*. 

In GhiringhcUi et al.^s paper,— the data were obtained 
for two different polarizations of the incoming photon: in 
the scattering plane and in the direction normal to the 
plane. We take the same geometry for our simulation 
with a plain MnOg octahedron. To first compute the 
multiplet, we consider an Mn ion with a formal neutral 
charge to avoid any shrinking of the wave-function and 6 
ions with a formal charge 2—, the Mn-0 distances 
are taken to be ~ 2.2 A and the octahedron is assumed 
to be undistorted. The results, scaled by 80% and with 
a constant broadening of 0.2 eV in the energy- loss di- 
rection, are shown in Fig [8] and compared to the actual 
experimental data. As shown in the lower panel of the 
figure, we have the possibility to investigate the polar- 
ization dependence for the outgoing photon this offers 
the opportunity to study the character of the different 
excitations depending on how each peak is affected by 
the selection rules. For instance, the peak close to an 
energy-loss of 4eV is more proeminent for crossed po- 
larizations. The effect of the polarization is also clearly 
visible on the elastic peak which weakens considerably : 
since the current effort in the study of transition metal 
oxides tends to investigate low-energy physics it would 
be of interest to take advantage of this effect to study 
the physics occurring in the vicinity of the elastic peak. 



V. CONCLUSION AND OUTLOOK 

In summary, we have described an easily accessible ap- 
proach for the computation of multiplet spectra arising 
in for soft X-ray spectroscopies of narrow band solids. 
The principal input specification is for the element where 
the core is to be excited and the core and valence con- 
figurations involved in the transition (for example: 2p^ 
3(P ) . The central field part of the present theory is done 
with a Dirac relativistic implementation for a simple den- 
sity functional, and yields thus the spin-orbit splitting 
from first principles. The electron-electron interaction 
amongst the open shell orbitals is obtained in the sub- 
space defined by the configuration specification. An in- 
tegral part of the approach is the treatment of arbitrary 
crystal fields. The specification of the crystal field is via 
straightforward specification of coordinates and charges 
for a small number of neighbor ions. The crystal field 
input remains easily oversecablc for arbitrary symmetry 
cases, without requiring expert knowledge in group no- 
tation. Further fundamental inputs arc the core hole 
lifetime, and polarization of the light with respect to 
the coordinate system used for the crystal field. With 
this minimal input, the method constitutes a first princi- 
ples approach, albeit with approximations. Comparison 
with experimental multiplet spectra is quite good usually 
without any fitting efforts. 



Nevertheless, the method can be tuned by a number of 
empirically adjustable parameters to fit experiment more 
accurately. The rationale for the empirical parameters 
are screening effects and coupling to energy bands that 
are not in the scope of the model. 

In its first version the method allows to calculate X-ray 
absorption spectra (XAS) and Resonant Inelastic X-ray 
Scattering spectra (RIXS). Comparison with experimen- 
tal spectra was presented in this article for some test 
cases. Planned extensions should include non ground 
state configurations in the RIXS. Extension to multi- 
plets in photo electron spectroscopy and circular mag- 
netic dichroism spectra is considered. Inclusion of charge 
transfer effects adjusted via empirical input parameters 
may also be considered in the near future. 

We plan to make this type of calculations widely ac- 
cessible in the near future via a web interface. 
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APPENDIX A: DOUBLE COUNTING 

The electron-electron interaction is given by the stan- 
dard Coulomb repulsion, with the following expansion: 

Ve-e = 1 r 

1=0 rn=-l ^ ^2 

with ri < r2 

(Al) 

Two terms appear in equation ([1]) : electrostatic 
interaction as well as eigenenergies coming form the 
solution of Dirac equation Cj. Through the self- 
consistentdetermination of the effective potential Veff 
appearing in the Dirac equation, these last terms contain 
already part of the electrostatic interactions: we thus 
should be careful and avoid double-counting. This can 
be achieved in subtracting the already counted part to 
Ea. (jAl[) by changing the sum X];^o ~* Si^i- The I = 



8 



term can be formally viewed as a ~ 1/r -term which is 
already treated in the potential of the Dirac equation. 
It is also worth mentioning that whenever a shell is 



full one can simply drop the electrostatic terms involving 
that shell since it only results in a common shift of the 
eigenenergies. 
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